Collective Charge Excitations below the Metal-to-Insulator Transition in BaVS,3 
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The charge response in the barium vanadium sulfide (BaVSs) single crystals is characterized by 
dc resistivity and low frequency dielectric spectroscopy. A broad relaxation mode in MHz range 
with huge dielectric constant ~ 10 6 emerges at the metal-to-insulator phase transition Tmi ~ 67 K, 
weakens with lowering temperature and eventually levels off below the magnetic transition T x ~ 
30 K. The mean relaxation time is thermally activated in a manner similar to the dc resistivity. 
These features are interpreted as signatures of the collective charge excitations characteristic for the 
orbital ordering that gradually develops below Tmi and stabilizes at long-range scale below T x . 

PACS numbers: 71.27. +a, 71.30.+h, 72.15.Nj, 77.22. Gm 



I. INTRODUCTION 

A variety of condensed matter systems with reduced 
dimensionality and strong Coulomb interactions between 
charges, spins and orbitals, display rich phase diagrams 
with novel forms of ordering phenomena 1 . These col- 
lective, broken symmetry phases like charge and spin- 
density waves, charge and orbital orderings have been 
the focus of intense scientific activity in recent years. 
Additionally, in transition metal compounds d-electrons 
experience competing forces: Coulomb repulsion tends 
to localize electrons, while hybridization with the ex- 
tended ligand valence states tends to delocalize them 2 -. 
The subtle balance makes these systems an excellent 
source for studying diverse metal-to-insulator (MI) phase 
transitions usually accompanied with drastic changes in 
charge, spin and orbital properties. In particular, the in- 
fluence of the orbital degeneracy and orbital ordering on 
electrical transport and magnetic properties is rapidly 
becoming a central issued. For example, in the classic 
transition metal oxide V2O3 X-ray anomalous scattering 
results have given a direct evidence that the orbital occu- 
pancy plays a central role in the physics of this system^. 
The spatial ordering of the occupancy of degenerate elec- 
tronic orbitals was shown to account for the anisotropic 
exchange integrals found in the antiferromagnetic insu- 
lator phase. Another outstanding issue is the question 
of collective excitations in the orbital ordered state. Al- 
though no theoretical model has been worked out until 
now, it is expected that the orbital degrees of freedom can 
either give rise to novel collective excitations, e.g. orbital 
waves, or can strongly renormalize other excitations*. 

The perovskite-type sulfide BaVS3 represents an ex- 
ceptional system to study the aforementioned phenom- 
ena and as such has attracted much attention in the last 
years. Face sharing VSg octahedra stack along the crys- 
tallographic c-axis and form the VS3 spin chains with 
one 3d electron per V 4+ site. The chains are sepa- 
rated by Ba atoms in the ab planes, yielding a quasi- 
one-dimensional (quasi-lD) structure. The unit cell at 



room temperature (RT) is a primitive hexagonal with 
two formula units. At 240 K the structure transforms 
into orthorhombic, but each chain keeps two equivalent 
V 4+ atoms per unit cell. The corresponding two electrons 
are shared essentially between two hybridized bands pro- 
duced by the crystal field splitting: a broad A\ g band 
derived from d z 2 orbitals overlapping along the c-axis, 
and a quasi-degenerate narrow E g \ band originating from 
e(i2g) orbitals with isotropic interactions via V-S-S-V 



bond: 



The filling of these bands is controlled by 



the on-site Coulomb repulsion U and local Hund's rule 
coupling J. The main effect is to bring the occupancies 
of A\ g and E g \ orbital closer to one another, which in the 
limit of strong correlations yields the occupancy for each 
of these orbitals close to half-filling. The spin degrees 
of freedom of the localized electrons and the coupling of 
conduction and localized electrons make the system ex- 
tremely complex. Their interplay results in the MI tran- 
sition at about 70 K and a magnetic transition at about 
30 K. 

In spite of a great deal of experimental efforts, no def- 
inite understanding has been reached yet on the detailed 
nature of the MI phase transition and the ground state in 
BaVS3 . Diffuse X-ray scattering experiments have shown 
that at Tmi a commensurate superstructure with the crit- 
ical wave vector q c = 0.5c* close to 2£^(Ai g ) sets in, 
preceded by a large fluctuation regime extending up to 
170 K 10 . This behavior is reminiscent of a Peierls transi- 
tion into a charge-density wave (CDW) state. However, 
the nature of the ground state is certainly more compli- 
cated since 2q c harmonic is found suggesting that the lo- 
calized e(t2g) electrons are also involved in the transition 
and order below T M i. Indeed, the related susceptibility, 
which follows the Curie- Weiss law from RT, exhibits an 
antiferromagnetic (AF)-like cusp at T MI and decreases 
belowii. There is no sign of magnetic long-range order 
down to T x « 30 K, where the incommensurate magnetic 
ordering is established 12 . The next unusual result at Tmi 
is a structural transformation from the orthorhombic to 
monoclinic with internal distortions of VSg octahedron 
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and tetramerization of V 4+ chains^ 3 .. The modulation 
of V distances in the superstructure can be described as 
a superposition of a 2kp and a 4/cf bond order waves 
(BOW), the latter component being primarily responsi- 
ble for the intensity change of the basic Bragg reflections. 
It is noteworthy that the tetramerization might be under- 
stood as an inherent feature of the Ikp Peierls transition, 
which reveals that a dimerization gap at high tempera- 
tures is negligible and that the broad A\ g band could be 
regarded as effectively one-quarter filled band. Finally 
and the most surprisingly, X-ray anomalous scattering at 
the vanadium if-edge revealed no charge disproportion- 
ation in the ground stated. A very intriguing interpreta- 
tion has been suggested involving the stabilization of two 
out-of-phase CDWs from d z 2 and efag) electrons, which 
implies an orbital ordering via an out-of-phase modula- 
tion of the occupancy of V sites by the d z 2 and efog) 
orbitals. 

In this paper, we address these important questions 
concerning the insulating state in BaVS3. Our results 
present a supporting evidence for an orbital ordered 
ground state and that the MI Peierls-like transition is 
also an orbital ordering transition. We find that a huge 
dielectric constant which arises in close vicinity of the 
MI transition dramatically decreases on cooling down to 
about 30 K and levels off below. We argue that such a 
behavior might be rather well explained in terms of an 
orbital ordering which sets in at the MI transition and 
develops the long-range order below the magnetic tran- 
sition. 



II. EXPERIMENTAL AND RESULTS 

dc resistivity was measured between RT and 10 K. 
In the frequency range 0.01 Hz-10 MHz the spectra of 
complex dielectric function were obtained from the com- 
plex conductance measured by two set-ups. At high 
frequencies (40 Hz-10 MHz) an Agilent 4294A preci- 
sion impedance analyser was used. At low frequencies 
(0.01 Hz-3 kHz) a set-up for measuring high-impedance 
samples was used^. The employed ac signal level of 
50 mV was well within the linear response regime. All 
measurements were done on single crystals along the 
crystallographic c-axis. The typical crystal dimensions 
were 3 x 0.25 x 0.25 mm 3 . 

An influence of extrinsic effects, especially those due 
to contact resistance and surface layer capacitance, was 
ruled out with scrutiny. In order to determine the qual- 
ity of contacts we have performed dc resistance measure- 
ments in the standard 4- and 2-contact configurations. 
Taking into account the difference in geometry between 
contacts in these two configurations, the 4-contact re- 
sistance can be scaled and subtracted from the 2-contact 
resistance in order to estimate contact resistance (R c ) vs. 
sample bulk resistance (R s )- Our analysis clearly shows 
that good quality contact samples whose dielectric re- 
sponse reflects the intrinsic features of the sample bulk 
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FIG. 1: Real and imaginary parts of the dielectric function 
of BaVS3 measured at three representative temperatures as a 
function of frequency with the ac electric field applied along 
the c-axis. The full lines are from fits by the generalized Debye 
expression (see Text). 



can easily be distinguished from bad samples by featuring 
the R c /R s ratio not larger than 5 in the metallic phase 
and in the vicinity of Tmi, and R c = 0.1i? s in the whole 
range of the insulating phase except close to Tmi- In con- 
trast, bad samples displayed R c /R s ratio of the order of 
100-1000 in the metallic phase and in the vicinity of Tmi, 
indicating large contact resistances, which concomitantly 
imply large contact capacitances throughout the whole 
measured temperature range. Consequently, the dielec- 
tric response registered in samples with bad quality con- 
tacts was a combination of the sample bulk and contact 
capacitance influence, the latter becoming dominant with 
increasing temperature. In this way, out of nine sam- 
ples studied, seven were discarded either due to bad con- 
tacts, or due to exceptionally low RRR at high pressure 
of 20 kbar. Two remaining good quality contact sam- 
ples with high RRR at high pressure have shown quali- 
tatively the same dielectric response, confirming that the 
observed response comes from the sample bulk. In this 
paper we present and discuss results obtained on one of 
these two single crystals. 

Fig. Q] shows frequency dependence of the complex di- 
electric response at three selected temperatures. A pro- 
nounced dielectric relaxation is observed. The screened 
loss peak e" centered at Tq 1 moves toward lower fre- 
quencies and smaller amplitudes with decreasing temper- 
ature. The main features of this relaxation are well de- 
scribed by the generalized Debye expression s(u>) — £hf = 
Ae/ [l + (iwto) 1 ""] , where Ae = eo — £hf (eo and chf 
are static and high-frequency dielectric constant, with the 
latter being negligible), tq is the mean relaxation time 
and 1 — a is the symmetric broadening of the relaxation 
time distribution function. 
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FIG. 2: Temperature dependence of the dielectric constant of 
the collective mode (upper panel), and dc resistivity and its 
logarithmic derivative (lower panel) in BaVS3. The arrows 
indicate the MI and magnetic transition temperature. 



Our results clearly demonstrate that a huge dielec- 
tric constant Ae is associated with the metal-to-insulator 
phase transition (Fig. [2]) . On decreasing temperature, a 
sharp growth of Ae starts in the close vicinity of Tmi and 
reaches the huge value of the order of 10 6 at Tmi = 67 K 
(Fig. O upper panel). This Tmi value corresponds per- 
fectly well to the phase transition temperature as deter- 
mined in the dc resistivity measurements, indicated by 
pronounced peaks at Tmi in the logarithmic derivative of 
resistivity (Fig. [2J lower panel). 

The observed dielectric relaxation would suggest a 
charge-density wave formation at Tvir^. The standard 
model of a deformable CDW pinned in a non-uniform 
impurity potential accounts for the existence of two 
modes, transverse and longitudinal 17 . The former cou- 
ples to the electromagnetic radiation and yields an un- 
screened pinned mode in the microwave region. Unfor- 
tunately, no microwave measurements on BaVS3 have 
been done yet. The longitudinal mode couples to an elec- 
trostatic potential and mixes in the transverse response 
due to the non-uniform pinning resulting in an over- 
damped low- frequency relaxation at Tq 1 due to screen- 
ing effects. The relaxation detected in our experiments 
bears two features as expected in the standard model. 
The first is that the relaxation time distribution is sym- 
metrically broadened, 1 — a ~ 0.8. The second is that 
the mean relaxation time To closely follows a thermally 
activated behavior similar to the dc resistivity tq(T) = 
T oexp(2A/2fc B T) oc p(T) (see Fig. 13]). r 00 w 1 ns de- 
scribes the microscopic relaxation time of the collective 
mode and the gap 2 A w 500 K, as found in the spec- 
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FIG. 3: The mean relaxation to of the collective charge mode 
(points) and the dc resistivity (full line) in BaVS3 as a func- 
tion of the inverse temperature. 



tra of the optical conductivity^. The dissipation can be 
naturally attributed to the screening due to single parti- 
cles originating from the wide A± g band. However, our 
discovery finds no consistent explanation within the stan- 
dard models, since the dielectric constant Ae displays a 
strong decrease below Tmi and levels off at temperatures 
below about 30 K, a behavior which significantly devi- 
ates from the one expected for a CDW condensate den- 
sity n: Ae(T) oc n(T). The decrease of Ae on moving 
between the MI transition Tmi and the magnetic tran- 
sition T x is substantial and amounts to two orders of 
magnitude. One possible explanation for this discrep- 
ancy is the very nature of the standard model for the re- 
sponse of the conventional CDW to applied electric fields 
in which the long-wavelength collective CDW excitation 
s.c. phason keeps the prominent role. We remind that 
this model is worked out for the incommensurate CDW 
in a random impurity potential, whereas the CDW in 
BaVS3 associated with the observed lattice modulation 
is commensurate with the order of commcnsurability, i.e. 
the ratio of the CDW and lattice periodicity, N = 4. 
However, the order of commcnsurability is not too high 
to impose the commensurability pinning and forbid the 
phason excitations— . Indeed, the experimental observa- 
tion of the broad relaxation, i.e. 1 — a « 0.8 indicates a 
randomness of the background structure. 

Another puzzling feature is brought by the dc electric- 
field dependent measurements. In the standard ID CDW 
compounds, the collective CDW mode couples to an ap- 
plied dc electric field and gives a novel contribution to 
the electrical conductivity-^. Our dc electric-field depen- 
dent measurements up to fields as high as 100 V/cm in 
the temperature range between 15 K and 40 K have only 
revealed a negligibly small non- linear conductivity, which 
emerges from the background noise. 
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III. DISCUSSION 

Our results indicate that the long-wavelength CDW 
collective excitations are frozen or strongly renormalized 
and that the collective excitations of different kind should 
be responsible for the observed dielectric relaxation. In 
the following we first address other possible causes and 
then offer the most plausible scenario for our results. 

First, we verify if the observed behavior of the di- 
electric constant might be due to the hopping conduc- 
tion which is known to arise in disordered systems with 
low dimensionality where Anderson localization yields 
the conductivity characterized by variable-range hopping 
(VRH) Mott law in dc limit, while in ac limit it fol- 
lows a power law dependence on frequency. Although 
BaVS3 may be regarded as a quasi-lD system, the hop- 
ping scenario does not seem realistic for the following 
reasons. First is that the frequency marking the on- 
set of the frequency-dependent transport is known to 
be roughly proportional to the dc conductivity, i.e. to 
the inverse of the dc resistivity 2 ^. The dc resistivity of 
BaVS3 at RT is about f CP 3 57cm and at lowest tempera- 
tures it is about 10 5 to f 6 Slcm. In the diverse systems 
with dc resistivities of similar orders of magnitude, the 
ac conductivity power law is observed only at frequen- 
cies above 1 MHz, whereas below I MHz an influence of 
hopping on dielectric dispersion is detected only for dc 
resistivities much higher than 10 10 Ocm i 21 ' 22 Indeed, in 
the case of BaVS3 a crude estimate for the crossover fre- 
quency yields values far above the frequency range where 
the dielectric response was observed 2 ^. The second result 
which excludes hopping comes from the observed optical 
spectral. Namely, a simple indication for existence of 
a hopping mechanism would be that the optical conduc- 
tivity is significantly enhanced in comparison with the 
dc conductivity, whereas in the case of BaVS3 the opti- 
cal conductivity, even at temperatures lower than Tmi, is 
at best comparable to the dc conductivity. 

Next possible origin of the observed dielectric response 
might be associated with a ferroelectric nature of the MI 
transition. Simple space group considerations indicate 
that below the MI transition the structure of BaVS3 is 
noncentrosymmetric with a polar axis in the reflection 
plane containing the VS3 chains. The symmetry of this 
low-tempcrature superstructure is Im, implying that the 
distortions of the two chains of the unit cell are out of 
phasei 3 -. Bond-valence sum (BVS) calculations of these 
X-ray data have indicated some charge disproportiona- 
tion at low temperatures. However, it was recently ar- 
gued by P. Foury-Leylekian^ 3 - that the BVS method over- 
estimated the charge disproportionation due to several 
reasons: a nonsymmetric V 4+ environment, thermal con- 
traction corrections which were not included and rather 
imprecise atomic coordinates which were used. Taking all 
this into account together with X-ray anomalous scat- 
tering result that showed negligible redistribution, not 
larger than 0.0 1 electron below Tmj^, we conclude that 
ferroelectricity cannot provide an explanation of the di- 



electric response in BaVS3. 

Finally, we address orbital ordering as a plausible 
ground state whose collective excitations might yield the 
observed dielectric relaxation. First we list arguments 
developed by Fagot et alJ^ who invoked an orbital or- 
dering associated with the MI phase transition in order 
to consistently explain structural data of BaVS3. Sup- 
porting results for an orbital ordering scenario at Tmi are 
an almost non-existant charge modulation in the insulat- 
ing phase, together with a qualitative structural anal- 
ysis of the VSg octahedron distortions, which reveals 
an out-of-phase modulation of the occupancy of V sites 
by the d z 2 and efog) orbitals. In particular, dominant 
E g i and A\ g occupancies are proposed for VI and V3 
sites respectively, while no definitive preferential occu- 
pancy was found for V2 and V4 sites. Furthermore, 
recent LDA+DMFT (local density approximation with 
dynamical mean-field theory) calculations performed in 
the monoclinic insulating phase of BaVS3 have qualita- 
tively confirmed an orbital-ordering scenario showing a 
V-site-dependent orbital occupancy and only minor, if 
any, charge disproportionation^. However, these calcu- 
lations suggest quite different orbital occupancies: it ap- 
pears that the (V3,V4) pair forms a correlated dimer with 
mixed A\ g and E g \ occupancy, while the VI and V2 ions 
bear major E g \ occupancy and negligible coupling. Fi- 
nally, the overall study indicates that although the local 
environment of the V site does not change substantially, 
the electronic structure turns out to be rather sensitive to 
change of temperature. The question arises what is the 
temperature dependence of the order parameter of or- 
bital ordering which starts to develop below Tmi and in 
what way it relates to the magnetic ordering. 51 V NMR 
and NQR measurements also suggested an orbital order- 
ing below Tmi that is fully developed only at T < Ty^i. 
The magnetic phase transition at T x is preceded by long- 
range dynamic AF correlations all the way up to Tmi 
and this phase bears features of a gapped spin-liquid-like 
phase. Mihaly et al<^ suggested that the lack of magnetic 
long-range order between Tmi and T x might be the con- 
sequence of the frustrated structure of a triangular array 
of V chains, which also prevents the orbital long-range 
order, so that the long-range spin and orbital orders can 
develop only well below Tmi- Finally, the AF static order 
below T x is not a conventional Neel phase: an AF do- 
main structure is suggested by the magnetic anisotropy 
measurements^. The existence of domains seems to be 
supported also by the muon spin rotation (/iSR) measure- 
ments, which showed an essentially random distribution 
of sizeable static electric fields below T x indicating an 
incommensurate or disordered magnetism 2 ^. 

Based on the considerations above we suggest the fol- 
lowing as the most plausible scenario. The primary order 
parameter for the MI phase transition is ID CDW in- 
stability and this CDW instability drives the orbital or- 
dering via structural changes involving a transformation 
from the orthorhombic to monoclinic with internal dis- 
tortions of VS6 octahedron and tetramerization of V 4+ 
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chains. The orbital order is coupled with the spin de- 
grees of freedom and drives the spin ordering into an 
AF-like ground state below 30 K. In other words, the or- 
bital ordering transition happens at Tmi, domains of or- 
bital order gradually develop in size with lowering tem- 
perature (concomitantly their number diminishes) and 
the long-range order eventually stabilizes below T x , al- 
beit domain structure persists. In this scenario we pro- 
pose that the collective excitations responsible for the 
observed features of the dielectric relaxation are short- 
wavelength ones, like charge domain walls^i in the ran- 
dom AF domain structure. Similar short-wavelength 
excitations associated with domain structure have pre- 
viously been invoked as the origin of dielectric relax- 
ation in diverse system s 28 ! 29 . The relaxation happens 
between different metastable states, which correspond to 
local changes of the spin configuration. The spin con- 
figuration is intimately connected with the charge and 
the orbital degrees of freedom as explained above. Since 
the dielectric constant is associated with the density of 
collective excitations, its anomalous temperature behav- 
ior below Tmi indicates that the relaxation-active num- 
ber of domain walls decreases with lowering temperature 
and eventually becomes well defined below T x . In other 
words, the dynamics of domain walls becomes progres- 
sively more restricted as the temperature lowers and it 
becomes constant below T x . 

IV. CONCLUSION 

In conclusion, we demonstrated the appearance of a 
huge dielectric constant associated with the metal-to- 



insulator phase transition in BaVS3 followed by a dra- 
matic decrease on cooling down to the magnetic tran- 
sition and leveling off below. We argue that the col- 
lective excitations whose dispersion we detect as broad 
screened relaxation modes arc not CDW phason exci- 
tations; rather they represent short-wavelength excita- 
tions of an orbital ordering, which sets in at the metal- 
to-insulator transition and develops the long-range order 
below the magnetic transition. Finally, BaVS3 repre- 
sents a beautiful example of the intricate interplay be- 
tween an orbital degeneracy on the one side, and spin 
and charge sector on the other side. This interplay needs 
to be taken into account in order to understand the ori- 
gin of the metal-to-insulator phase transition and low- 
temperature phases in the transition metal compounds 
in general. Further work on the theoretical and experi- 
mental fronts is needed to demonstrate directly the ex- 
istence of an orbital order in BaVS3 and the associated 
superstructure. 
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